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ABSTRACT: The swine leukocyte antigen class II molecules are possibly associated with the induction of protective immunity. The 
study described here was to investigate the relationship between polymorphisms in exon 2 of the swine DQA gene and piglet diarrhea. 
This study was carried out on 425 suckling piglets from three purebred pig strains (Large White, Landrace and Duroc). The genetic 
diversity of exon 2 in swine DQA was detected by PCR-SSCP and sequencing analysis, eight unique SSCP patterns (AB, BB, BC, CC, 
CD, BD, BE and DD) representing five specific allele (A to E) sequences were detected. Sequence analysis revealed 21 nucleotide 
variable sites and resulting in 12 amino acid substitutions in the populations. A moderate level polymorphism and significant deviations 
from Hardy- Weinberg equilibrium of the genotypes distribution were observed in the populations (p<0.01). The association analysis 
indicated that there was a statistically significant difference in the score of piglet diarrhea between different genotypes, individuals with 
genotype CC showed a lower diarrhea score than genotypes AB (0.98±0.09), BB (0.85±0.77) and BC (1.25±0.23) (p<0.05), and 
significantly low than genotype BE (1.19±0.19) (p<0.01), CC genotype may be a most resistance genotype for piglet diarrhea. (Key 
Words: SLA, DQA Gene, Genetic Polymorphism, Piglet Diarrhea) 



INTRODUCTION 

The Large White, Landrace, and Duroc are three 
important commercial pig breeds in the intensive breeding 
pig farms and are well-known for their characteristics of 
high prolificacy, fast growth and high lean rates. Due to 
large scale intensive production systems, pigs are exposed 
on various pathogens that cause infectious diseases and 
pathological symptoms (Lunney et al., 2009). 

Piglet diarrhea is one of the serious threats to economic 
success in the swine industry, which causes increased 
mortality rates, poor growth and additional medical costs. It 
is estimated that 50% of piglet mortality is due to diarrhea 
deaths (Morris et al., 2002). Generally, the incidence of this 
disease is controlled by the combined function of genetic 
factors and pathogens such as contagious and 
environmental pathogens and poor management practices. 
An effective stimulation of the adaptive immune system is 
based on the activation of T lymphocytes via their 
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recognition of pathogen-derived peptides presented by 
major histocompatibility complex (MHC) molecules to the 
respective T-cell receptors. The swine MHC, which codes 
for a series of extremely polymorphic cell-surface 
glycoproteins, is referred to swine leukocyte antigens (SLA) 
which function is to capture antigen peptides as well as 
present them to T-lymphocytes, thereby influencing host 
challenge to infectious disease (Lunney et al., 2009; Patch 
et al., 2011). 

The associations of porcine SLA complex with swine 
immune responses, as well as disease resistance and 
susceptibility, have led to a number of comprehensive 
studies. SLA alleles and haplotypes have been associated 
with difference in swine immune response to various 
infectious agents (Mallard et al., 1989; Renard and Vaiman, 
1989; Shinkai et al., 2012). In the NIH miniature pigs, SLA 
haplotypes have been shown to affect antibody responses to 
sheep erythrocytes and cell-mediated immune responses to 
tuberculin-purified protein (Mallard et al., 1989). Tissot et 
al. (1989) also found a relationship of haplotype B with a 
higher susceptibility of swine cutaneous malignant 
melanoma in Sinclair miniature swine. A possible SLA 
influence on the inherent disease susceptibility and 
resistance in PRRS has also been reported (Lewis et al., 
2007; Molina et al., 2008). Yang et al. (2012) reported that 
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polymorphism of the SLA-DRA gene has significant 
influence on diarrhea of piglets. All of those works 
indicated that a detailed analysis of SLA genes is essential 
to advance our understanding of the processes implicated in 
immune responses and the effects of intensive selection, 
such as via strict breeding programs. 

Generally, the extensive polymorphism of MHC genes 
is expected to broaden the array of antigenic peptides that 
available for T-cell immunological recognition in 
stimulating and regulating immune surveillance. This 
extremely polymorphism, is essentially concentrated in the 
first domain (exon 2), in a region which forms the floor and 
wall of the peptide-binding groove of the class II genes 
(Smith et al., 2005). 

The SLA class II genes expressed on the surface of 
professional antigen presenting cells are noted as SLA-DRA, 
SLA-DRB1, SLA-DQA, and SLA-DQB1. SLA-DRB1 and 
SLA-DQB1 loci display a very high degree of 
polymorphism, although the SLA-DQA gene has been 
reported as a moderately polymorphic gene, and SLA-DRA 
was reported to exhibit only limited polymorphism (Lunney 
et al., 2009; Ho et al., 2009; Ho et al., 2010). Research 
utilizing SLA class II genes primarily focus on exon 2 of 
the DRB and DQB because these loci code for the peptide 
binding region (PBR) which is in direct contact with the 
processed peptide (Babik et al., 2005; Luetkemeier et al., 
2009; Park et al., 2010). However, there is currently no 
study on the relationship between SLA-DQA gene and piglet 
diarrhea. The swine DQA gene encodes the a chain of DQ 
antigen which is closely linked to DRB. The full-length of 
the genome sequence of this gene contains 5.5 kb that 
encodes for a polypeptide of 255 amino acid residues, 
which include four exons, with exon 1 encoding signal 
peptide region, with exon 2 and 3 encoding the 
corresponding extracellular al and a2 domains, and exon 4 
encoding both transmembrane and cytoplasmic domains 
(Lunney et al., 2009). 

Consequently, in the present study we characterized the 
genetic variation of the exon 2 of SLA II DQA locus from 
Large White, Landrace and Duroc pig strains and 
investigate the correlation between genotypes and piglet 
diarrhea, to identify the suitable candidate markers 
correlated with animal resistance potentially for use in 
genetic strategies to breed disease resistance pigs. 

MATERIALS AND METHODS 

Samples and data collection 

The animals used in this experiment were from Linze 
Xinghua Pig Breeding Farm (Zhangye, Gansu, China). A 
total of 425 pure piglets were examined, including Landrace 
(n = 146), Large White (n = 148) and Duroc (n = 131) pigs. 

Experimental piglets were born in June 2010 and the 



experiment conducted under forage and feeding 
management conditions and disease control in accordance 
with the Guide for Chinese Feeding Standard of Swine and 
approved by the institute of Ministry of Agriculture Feed 
Industry Center, China Agricultural University (Beijing, 
China). All piglets were monitored two times a day during 
the entire suckling period (from birth to weaning, 0 to 28 d) 
and assigned a daily score based on visual analysis of 
symptoms traits as follows: 0- normal, solid feces; 1- slight 
diarrhea, soft and loose feces; 2- moderate diarrhea, semi- 
liquid feces; 3- severe diarrhea, liquid and unformed feces 
(Kelly et al., 1990). Genomic DNA was extracted from 
swine ear tissues samples by standard phenol-chloroform 
extraction procedure (Sambrook and Russell, 2001) and 
quantified by using a Nanodrop ND-1000 
spectrophotometer (Thermo Fisher Scientific, USA). 

Genotyping by polymerase chain reaction-single-strand 
conformation polymorphism (PCR-SSCP) 

PCR was performed to amplify a 378 bp fragment of the 
SLA-DQA gene which encompassed 81 bp intron 1, 249 bp 
exon 2 and 48 bp intron 2 using forward primer SLA- 
DQAe2F (5 -AGTCAAGTTCTCTTGTCACT-3') and SLA- 
DQAe2R (5 ' -TGTG A ACGGGTAG ATTCTGT-3 ') . The 
primers were designed according to the published SLA- 
DQA sequence in GenBank (accession No. AY303988). 
PCR analyses were carried out in total of a 25 |il volume 
with 100 ng genomic DNA, 10 pmol of each primer, 2.5 
mM dNTPs, 5 U of Taq DNA polymerase (TaKaRa, Dalian, 
China), and lOxreaction buffer (including 15 mM Mg 2+ ). 
PCR conditions included an initial denaturing of 2 min at 
94°C, followed by 30 cycles of a three-step process of 
denaturation at 94 °C for 20 s, 55 °C for 30 s and extension 
at 72°C for 30 s, followed by a final extension at 72°C for 
10 min. PCR products were detected on a 2% agarose gel 
electrophoresis. 

All samples were screened using SSCP analysis to scan 
for polymorphism in the amplified region. About 2 |il PCR 
product was mixed with 6 |il of the denaturing solution 
(98% formamide, 0.025% bromophenol blue, 0.025% 
xylene cyanole, and 10 mM EDTA), denatured for 10 min 
at 98°C, and immediately chilled on ice for 10 min. 
Denatured DNA samples were placed on 17 cmxl7 cm, 
10% acrylamide: bisacrylamide gel (29: 1) gels. 
Electrophoresis was run in 1 xTBE buffer of 30 min (250V) 
at 4°C and then 15 h (180V) at room temperture. Gels were 
silver-stained according to the method of Byun et al. (2009). 

Cloning and sequencing 

To determine the sequence variation in the swine DQA 
gene, PCR products of individuals that showed different 
SSCP patterns were purified using DNA Fragment Quick 
Recover Kit (TaKaRa, Dalian, China). The purified PCR 
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products were ligated into pMD®19-T vectors (TaKaRa, 
Dalian, China) and transformed into the Escherichia coli 
DH 5 a strain. Recombinant clones were detected by 
blue/white screening, between 10 to 15 insert-positive 
clones for each transformation were picked and incubated 
overnight in terrific Luria-Bertani broth at 37°C, in a 
vibrating rotary incubator (200 rpm). Clones were screened 
using the PCR-SSCP approach as described by Zhou et al. 
(2008), and only those clones for which the SSCP patterns 
matched those of the corresponding genomic DNA were 
selected and subsequently sequenced by Sangon Biotech 
(Shanghai, China) Co., Ltd. 

Statistical analysis 

Sequence alignment, translation and comparison were 
performed using MEGA5 (Tamura et al., 2011); genotype 
and allele frequencies, observed heterozygosity {He), and 
testing for deviation from Hardy-Weinberg equilibrium 
were calculated by the POPGENE version 1.32 (1997); 
polymorphism information content (PIC) was calculated by 
Botstein's methods (1980). 

The average rates of nonsynonymous (d N ) and 
synonymous (d s ) substitution for antigen-binding sites 
(ABS) and non-antigen binding sites (non-ABS) were 
calculated according to Nei and Gojobori (1986) method 
with Jukes-Cantor (1969) correction, using the package 
MEGA5 (Tamura et al., 2011). The porcine ABS was 
defined by assuming congruence with human MHC sites 
(Brown et al., 1993). 

The effects of genotype on piglet diarrhea were 
estimated using the general linear model (GLM) procedure 
of SAS software (SAS, 1999) according to the following 
statistical model: 

Y ijkim = H+Bi+Sj+Gii+eyH 

Where is the observation value of piglet diarrhea 
score, |i is the overall population means, Bj is the fixed 
effect of breed, Sj is fixed effect of gender, G k is the effect 
of genotype, and e is the random residual effect. 
Preliminary analysis also included the fixed interaction 
effects of breed, sex and genotype, however, these 
interaction effects were subsequently removed because it 
did not have a significant effect. 

RESULTS 

Polymorphism analysis in exon 2 of the SLA-DQA gene 

A 378 bp amplified fragment of the SLA-DQA gene was 
obtained. Based on the PCR products, SSCP was used for 
initial screening of samples, a total of 8 unique SSCP 
patterns (AB, BB, BC, CC, CD, BD, BE and DD) were 
detected among 425 pigs (Figure 1). Genotypes AB, BB, 



AB BC BC BE BB DD CD BD BD CC 




Figure 1. The genotyping results of the SLA-DQA gene exon 2 by 
PCR-SSCP. 



BD, BE and BC were detected in all three pig strains, 
genotype CC was only detected in the Landrace breed and 
genotypes CD and DD were only detected in Duroc breed. 
Homozygous individuals were found for alleles B, C and D, 
but alleles A and E only appeared in heterozygous 
individuals. Allele B was the most common allele, with 
frequencies of 0.61, 077, and 0.69 in Large White, Landrace 
and Duroc pigs, respectively. 

Significant deviations from Hardy-Weinberg 
equilibrium were detected in the genotypes distribution in 
Large White, Landrace and Duroc pig populations (p<0.01). 
The observed levels of heterozygosity ranged from 0.39 to 
0.54, and the PIC for the studied populations ranged from 
0.37 to 0.47 (Table 1). According to the protocol of Vaiman 
et al. (1994), this locus is moderately polymorphic in the 
three populations (0.25<P/C<0.50). 

Table 1. Genotype and allele frequency, population genetic 
parameters of the SLA-DQA exon 2 in three pig strains 
Breed Large White Landrace Duroc 

n 146 148 131 



Genotype frequency 



AB 


0.60 (87) 


0.14 (20) 


0.16(21) 


BB 


0.22 (32) 


0.60 (89) 


0.46 (60) 


BD 


0.09 (13) 


0.04(6) 


0.04 (5) 


BE 


0.05 (8) 


0.10(15) 


0.14(18) 


BC 


0.04 (6) 


0.07 (10) 


0.14(18) 


CC 


0.00 (0) 


0.06 (8) 


0.00 (0) 


CD 


0.00 (0) 


0.00 (0) 


0.05 (6) 


DD 


0.00 (0) 


0.00 (0) 


0.02 (3) 



Allele frequency 



A 


0.30 


0.07 


0.08 


B 


0.61 


0.77 


0.69 


C 


0.02 


0.09 


0.09 


D 


0.04 


0.02 


0.06 


E 


0.03 


0.05 


0.07 


x 2 (HWE) 


59.29* 


57.37* 


43.74* 


He 


0.54 


0.39 


0.49 


PIC 


0.47 


0.37 


0.46 



x 2 (HWE) = Chi-square test for Hardy-Weinberg equilibrium. 

He = Observed heterozygosity. PIC = Polymorphism information content. 

* Hardy-Weinberg disequilibrium, p<0.0 1 . 
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3898 3942 396T 

AY3039S8 CCGACCATGTTGCCTCCTAIGGCTTAAATGTCTACCAGTCTTACGGTCCCAGCGGCTATTATACCCATGAATTTGATGGCGACGAGGAAT 

SLA-D0A'E2s 

SLA-D6A*E2b 

SLA-D0A*E2c m > T c • • • 

SLA-DQA'- 2d * C... 

SLA-DQA*E2e- 



AY3039S8 TCTATGTGGACCTGGAGAAGAAGGAGACTGTCTGGCAGCTGCCTCTGTTTAGCAAATTTACAAGITTTGACCCGCAGGGTGCACTGAGGA 

SLA-DQA*E2a 

SU-D0A*E2b C C TG 

SU-D0A*E2c" C A 

SU-WMBd C C TG 

SU-DdA*E2e' C.A C TG 

4078 4122 



AY303998S ACATAGCTACGTTAAAACATAACTTGAACATTGTGACTAAACGCTCCAACAACACCGCGGCTGTCAATC 

SLA-DOA*E2a A 

SLA-DOA*E2b - T c 

SLA-DQA*E2c* GT.GG T..C CC A A 

SLA-DQA*E2d * 

SLA-DQA*E2? -T C 



Figure 2. Alignment of the SLA-DQA exon 2 allele nucleotide sequences refer to the sequence AY303988. Sequences correspond to 
positions 3898-4146 were the complete exon 2. Identical nucleotides are replaced by • and the substitutions are written as letters. Newly 
identified sequences are marked with * which have assigned GenBank accession numbers. 



+ + + + ++ + +++ ++ + ++ +■ ++ + 

AY303988 DHVASYGLNVYfiSYGPSGYYTHEFDGDEEFYVDLEEKETVWgLPLFSKFTSFDPQGALRNIATLKHNLNIVTKRSNNTAAVN 

SlA-DQA'Uu 

SLiDQA 'E2b p E 

SLl-DQA'E2c 3 Q 3 VG.LT...L 

SU-DQA"E2d R S P E 

SLA-DQA'E2t E P E 

Figure 3. Alignment of the SLA-DQA exon 2 allele amino acid sequences refer to the sequence AY30398. Putative antigen-binding sites 
(ABS) are marked with plus which have been defined according to Brown et al. (1993). 



The nucleotide characterizations of these 
polymorphisms were revealed by cloning and sequencing. 
Comparison with the existing alleles in the Immuno 
Polymorphism Database (IPD)-MHC SLA sequence 
database (http://www.ebi.ac.uk/ipd/mhc/sla/) found two 
new mutations (c.3994T>A and c.4065G>A) for the first 
time in swine, which formed two new alleles. The new 
allele sequences have been deposited in GenBank database 
(http://www.ncbi.nlm.nih.gov/nuccore/) under accession 
number JX840404, JX840405. 

Alignments of nucleotide and amino acid sequences 
showed that twenty-one of 249 (8.43%) nucleotide (Figure 
2) and twelve of 82 (14.63%) amino acid variation sites 
(Figure 3) were identified in exon 2 of the swine DQA gene. 
We estimated the amino acid variation separately for PBR 
and non-PBR, the nonsynonymous substitutions detected 
cause changes at 12 amino acid sites, 6 (31.57%) were 
included in the 19 PBR (marked with+in Figure 3) sites, 
whereas only 6 (9.52%) of the 63 amino acid positions were 
polymorphic in the non-PBR. The rate of synonymous 
substitutions among DQA alleles occurred at higher 
frequencies than nonsynonymous substitutions at all sites as 
well as for the ABS and non-ABS in this study (d N /d s <1, 
purifying selection) (Table 2). 



Effects of genotypes on piglet diarrhea 

Statistical results showed that a significantly difference 
at swine DQA exon 2 locus was observed between genotype 
effect and piglet diarrhea (p<0.05), whereas breed and sex 
showed little influence on piglet diarrhea (p>0.05). The 
least square means along with standard errors for piglet 
diarrhea score among different genotypes are given in Table 
3. The genotypic linear contrasts suggested that diarrhea 
score of genotypes CC (0.36+0.19) was extremely lower 
than that of genotype BE (1.19+0.19) (p<0.01), and 



Table 2. The average rates of synonymous substitutions (dS) and 
nonsynonymous substitutions (dN) among 5 SLA-DQA exon 2 
alleles 



Positions 


N 


Nv 


dN 


dS 


dN/dS 


Combined 


82 


12 


0.034 


0.055 


0.618 








(0.008) 


(0.009) 




PBR 


19 


6 


0.085 


0.109 


0.800 








(0.028) 


(0.023) 




Non-PBR 


63 


6 


0.019 


0.041 


0.463 








(0.003) 


(0.010) 





N = Numbers of encoded amino acid. 
Nv = Numbers of variable amino acid. 

Standard errors (in parentheses) were obtained through 1,000 bootstrap 
replicates. 
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Table 3. Association of different genotypes of the exon 2 of SLA- 
DQA gene with piglet diarrhea score 



Genotype 


Samples 
number 


Least square means±standard errors 


AB 


130 


0.98±0.09 a 


BB 


173 


0.8510.77" 


BD 


43 


0.80±0.15 ac 


BE 


26 


1.19±0.19 A 


BC 


18 


1.25±0.23 a 


CC 


27 


0.36±0.19 Bbc 


CD 


5 


1.12±0.43 ac 


DD 


3 


1.23±0.56 ac 



Values with different superscripts within the same line differ 
significantly at p<0.01. 

b ' c Values with different superscripts within the same line differ 
significantly at p<0.05. 



significantly lower than genotypes AB (0.98+0.09), BB 
(0.85±0.77) and BC (1.25+0.23) (p<0.05) (Table 3). The 
animals with genotype DD were excluded when association 
analysis was carried out because the sample size was less 
than three. 

DISCUSSION 

Genetic diversity of the SLA-DQA gene 

One of the most remarkable features of the MHC genes 
is the extremely high degree of genetic polymorphism 
within loci. Extensive genetic studies have been carried out 
on the MHC II DQA gene of various species (Hickford et al., 
2007; Janova et al., 2009); 24 previously DQA alleles and 
13 new alleles were identified in the genus Equus by Janova 
et al. (2009), and 37 ovine DQA haplotypes were detected 
in 520 sheep by Hickford et al. (2007). Variation in swine 
the DQA gene is well characterized, 20 alleles have been 
identified with a moderate degree of polymorphism 
(Lunney et al., 2009; Ho et al., 2010). In this study, we have 
successfully characterized the coding sequence of swine 
DQA exon 2, a total of twenty-one nucleotide variation sites 
were identified, defining 5 alleles. In particular, with the 
contribution of two novel alleles was identified. These 
indicate that the SLA-DQA gene, which is similar to other 
vertebrates, is one of the most polymorphic genetic loci, 
with high allelic diversity and high levels of sequence 
divergence among alleles. 

The distribution of all genotypes had significantly 
deviations from Hardy- Weinberg equilibrium (p>0.05) in 
Large White, Landrace and Duroc pig populations. Five 
heterozygous genotypes (AB, BC, BE, CD and BD) and 
three homozygous genotypes (BB, CC and DD) were 
present in this study. This may attributable to the difference 
of survival in an environment with selection pressure from 
antigens and increased other environmental influences 



acting to dilute the impact of genetic variation (Andersson 
et al., 2012). The advantage for heterozygous individuals 
may be due to their wider antigen recognition spectrum than 
homozygous against a wide variety of pathogens (Sommer 
et al., 2005), while missing homozygous genotypes may 
have been eliminated in the breeding process (Beuzen et al., 
2000). 

He and PIC are considered to be suitable parameters to 
estimate the genetic variation of populations. The average 
diversity in SLA II genes has been shown to have a high 
polymorphism with He (0.628) and PIC (0.581) (Yu et al., 
2012). The He and PIC of our study were greater than 0.37 
in all three populations, this concurs with the results of 
Lunney et al. (2009) that indicate SLA-DQA diplays a 
middle level of polymorphism, compared to previous swine 
SLA studies by Yu et al. (2012). 

A popular explanation for the extensive polymorphic 
nature of MHC genes is that it is driven and maintained by 
balancing selection. The classical way to investigate 
whether an MHC gene has been subject to balancing 
selection is to investigate the ratio between non- 
synonymous and synonymous substitutions (d N /d s ) in the 
PBR. Intriguingly, in the present investigation a low 
proportion of d N and d s at PBR suggests that this gene 
actually has been subject to purifying selection. This is in 
contrast to the previous DQA studies in other vertebrates 
species (Zhu et al., 2007; Goiiy de Bellocq et al., 2009; 
Kamathand and Getz, 2011). Through purifying selection, 
non-synonymous substitutions that probably affect normal 
gene function will usually be eliminated (Hughes and 
Yeager, 1998). 

Association between genotypes and piglet diarrhea 

Previously studies have established that DQA gene is 
implicated in human atopic dermatitis, Puumala virus and 
Cowpox virus infections in bank voles (Deter et al., 2008; 
Kiyohara et al., 2008). In pigs, SLA alleles and haplotypes 
had been confirmed to be associated with some immune- 
related diseases (Tissot et al., 1989; Chen et al., 2005). 
Renard and Vaiman (1989) indicated that mortality rate 
before weaning caused by diarrhea was linked to SLA 
haplotypes. 

The SLA class II DQA gene encodes cell-surface 
glycoproteins that preferentially recognize foreign antigens 
and, subsequently, alter disease and vaccine responsiveness 
(Smith et al., 2005). However, there is currently no study on 
the relationship between SLA-DQA gene and piglet diarrhea, 
however, a study by Yang et al. (2012) reported the diarrhea 
score of SLA class II DRA genotypes AA (1.10+0.19) and 
BB (1.05+0.14) significantly higher than AC (0.32+0.25) 
and CC (0.06+0.12). In the present study, we also observed 
a statistically significant difference in the score of piglet 
diarrhea between different genotypes, genotype CC 



1070 



Yang et al. (2013) Asian Australas. J. Anim. Sci. 26:1065-1071 



(0.36+0.19) indicated a lower diarrhea score than genotypes 
AB (0.98±0.09), BB (0.85+0.77), BC (1.25+0.23), and BE 
(1.19+0.19). This finding supports the previous 
observations by Yang et al. (2012) and Renard and Vaiman 
(1989). 

Amino acid mutations in the coding region have a 
significant impact on the stability of proteins, which largely 
depend on the chemical and physical properties of the 
mutated amino acid, these include mutation position 
(whether in functional domain) and charge properties 
(Wang and Moult, 2001). Six of the 12 polymorphic amino 
acid sites were included in the peptide binding region in our 
study, which might play a crucial role in immunological 
processes modulating disease pathogenesis. Especially, in 
the case of c.4086A>G and c.4087C>T, as linked AC G 
nucleotide variation into GT G of allele C, changes the 
encoded amino acid from threonine to valine; and 
c.4089A>G and c.4090C>T, as linked TTA nucleotide into 
GG A of allele C, changes the encoded amino acid from 
leucine to glycine. These are notable as both caused 
changes of amino acid properties from hydrophobic and 
polar to lyophobic and nonpolar and therefore may affect 
protein function or stability causing variation in genotype. 
Our statistical analysis revealed CC genotype showed a 
lower diarrhea score than other genotypes, indicating that 
those four nucleotide mutations are possible functional 
mutations related to piglet diarrhea and genotype CC may 
be a most resistance genotype for piglet diarrhea. 

To our knowledge, this is the first report describing a 
significant association of the variation in swine DQA with 
piglet diarrhea. The results reported here provide a 
theoretical foundation for further pig anti -disease resistance 
breeding research. Considering piglet diarrhea is influenced 
by multiple factors, such as parasite disease biology, genetic 
background, and sample size and especially MHC typing 
method, further works (sample amplification and DNA 
typing) should be carried out to facilitate future 
understanding the associations with specific alleles and 
SNPs. This will provide more meaningful information with 
possible application in future pig disease resistance 
breeding programs. 
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